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the simplicity and ease of application of the theory, it is a valuable aid in 
soils work. So that its use may become more widespread, description of the 
theory and its use is illustrated in the examples that follow. 

Professor Pavlovsky in his theoretical development considered only em- 
bankments with impervious foundations, because from the hydraulic point of 
view an embankment on an impervious foundation giv2s the highest elevation 
of the phreatic line. However, the theory is also app.oximately applicable to 
dams on pervious foundations. In the latter case, the seepage water appears 
in the foundation as well as in the fill itself. The rate of percolation and the 
volume of flow through the compacted fill will differ from the rate and vol- 
ume through the foundation largely because of the difference in type and prop- 
erties of the material. A dam designed for an impervious base, but to be 
constructed on a porous one, will have a line of saturation lower than assumed, 
and the seepage pressure on the downstream slope will be less. This will add 
stability and safety to the embankment section or will permit economic modi- 
fications in the design. The derivation of the formulas for the seepage flow 
through the homogenous dam and for locating the theorectical line of satura- 
tion was simplified by considering a section of an earth fill of unit thickness, 
divided into three characteristic parts which are shown in Figure 1 and are 
called upstream, central, and downstream parts, respectively. The general 
hydraulic properties of the percolation water through each of these parts are 
analyzed separately. 

In the upstream part of the earth dam, after water penetrates the upstream 
slope, the water starts to flow into this part of earth fill (as shown in Figure 1, 
between the upstream slope and Section W-W). During percolation of water 
through this part, the water meets resistance which requires some energy to 
overcome. Therefore, in the second part of its path, in the central part, it 
has much less energy. As in the motion of ground water, the kinetic energy 
is negligible compared with the potential energy. The potential head referred 
to datum 0-0 is represented by Hy, the depth of the water on the upstream 
side. Therefore, the loss of head a (in the upstream part) is expressed by 
the difference between H,, and h, where h is the remainder of the total poten- 
tial energy of the water after percolating through the upstream part of the 
earth dam. 

From Section W-W (Figure 1), the percolating water does not have a uni- 
form moticn through the soil, and the line of saturation in the central part of 
the dam will be a descending curve and will intersect the surface of down- 
stream slope at point Mp. 

From hydraulic considerations, the position of this point will be located 
above the tail water surface. The position of M, is utilized as a starting 
point in analyzing the stability of the downstream slope of an earth fill as 
will be seen in the second part of this article. 

In studying the general hydraulic properties of the downstream part of 
dam, the following statements apply: The potential energy of water for the 
entire dam is Hy - ho. If the outlet point, M5, of the line of saturation is at 
the elevation of the tail water (Figure 1), the potential energy for the down- 
stream part is zero, and there would be no flow in the downstream part. But 
»ince there is flow in the downstream part, it must be concluded that the out- 
let point M, must be above the level of the tail water. Hence a, is introduced 
to represent the height of the outlet point M, above the elevation of the tail 
water, and the following equation can be written 
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From the above statement, the following conclusion can be drawn: The down- 
stream part of the earth dam has two zones of percolation flow: that is, the 
upper zone, with a discharge of percolation flow q,, which goes to the atmos- 
phere, and the second, below the first zone under level of the tail water with 
a discharge q9. Therefore, the total discharge through the section of the 
earth dam of unit thickness is 

= 4; + 

In previous discussion, the case was considered in which the depth of tail 
water h, >0, but in practice h, is usually zero; that is, there is usually no 
water at the downstream side of the earth dam, In this case, the mathemati- 
cal treatment concerning the upstream part, and the central part of the earth 
dam will be the same as in previous case. But in the downstream part, the 
position of the outlet point M, of the line of saturation will change and will 
be at an elevation a, above the base line (Figure 2). This elevation is neces- 
sary to overcome the hydraulic resistance during the percolation of water 
through the downstream part of the dam. 

Concerning the general hydraulic properties of the percolation through an 
earth dam, it is necessary to consider the direction of the flow lines at the 
entrance to the body of the dam at the upstream part. The upstream slope 
represents the equipotential line. It is known that the equipotential lines are 
normal to the flow lines. So it is evident that the flow lines in an earth dam 
must be normal to the surface of the upstream slope at the time of entrance 
to the upstream part. During the percolation flow in the upstream part, the 
flow lines will change their direction and at the time of entrance to the bound- 
ary of the central part, they will generally choose a nearly horizontal direc- 
tion. This conclusion applies to the line of saturation itself, which represents 
the upper flow line of the percolating water. From description of the main 
scheme of percolation of water through a homogenous earth dam on an imper- 
vious foundation, the flow in the body of the dam can be considered to be a 
hydraulic system with three divisions; i.e., the upstream part plus the central 
part plus the downstream part. These three parts represent one hydraulic 
unit, and through each part the same amount of water, q, percolates. Deriva- 
tions of the final usable equations pertaining to these problems are presented 
in the appendix of this paper. However, the study and more detail concerning 
the other types of earth fill section may be found in “Technical Memorandum 
No. 383,” Bureau of Reclamation. 

Applying the Dupuit formula and Darcy’s law for the flow of water through 
a homogenous earth dam, the general equation for the line of saturation for 
the central part of the dam was derived, the mathematical expression of 
which is: 


=n? - (1) 


where, 


y = vertical ordinate 

X = positive distance from W-W to point considered 

k = coefficient of permeability 

q = percolation discharge per unit length of the dam and 
h = the elevation of point Mj), above the base line 


For the percolation discharge through the body of an earth dam, three 
equations were developed by Pavlovsky; i.e., for the upstream part, the 
central part, and the downstream part. When h,> 0 as in Figure 1 these 
equations have the following expressions: 
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These three equations include four unknowns: 
q, h, and S. 


However, one of the four unknowns, S, can be easily determined from the 
hydraulic boundary of the earth dam (Figure 1) and 


S=b+m, [Hy - @, +h] 5) 


Equations (2) through (5) determine the ordinates of the main points of the 
line of saturation, M, M,, and M, (Figure 1), and together with Equation (1) 


the intermediate points between M, and M,: 


From study of equations (2 to 4) the position of the line of saturation in an 
earth dam does not depend upon the coefficient of permeability of the earth 
material from which the dam is constructed, but depends only on the hydraulic 
boundaries of the structure. 

After simultaneous solution of the above equations, determining the value 
k, and knowing the coefficient of permeability, k, the percolation discharge 
per unit length of dam can be computed. 

The first part of the line of saturation in the upstream part of the dam, 
between the point M and M,, cannot be plotted by points, but the above meth- 
od gives the position of the two points, M and M,, where the first point is 
determined by the intersection of the upstream surface water with the slope 
of the dam and the second by the value of h at M. 


The line of saturation between the points M, and M, can be plotted com- 
pletely using Equation 1, and the curve M to M, must give a smooth connec- 
tion to have a common tangent at point M, with the curve from M, to M,: 


The tangent at point M, of the line of saturation forms an acute angle (6) 


1 
with the axis of the abscissa and is described by 


tan@ =q_ 
kh (6) 

This explanation is sufficient for practical plotting of the line of saturation 
in the body of the homogenous dam because a more accurate plotting of the 
line of saturation is not necessary in designing the section of an embankment. 

For simplification and ease of application of the percolation flow equations, 
the following designation is introduced in Equation 4. 


A= 1 + log, a 


(7) 


= 1 
a= 35 ho)” | (3) 
= 
(4) 
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Then Equation 4 will be 


A 


Substituting the above expression in Equation 3 and solving for h we have 


[ 28a 2 
— A+ (a, (9) 
1 a 


and substituting the expression for S from Equation 5, Equation 9 becomes 
b 
h = 2a,A + Ha (a, + h,)| (10) 
Equation 2 can be simplified by having 
(11) 


then Equation 2 becomes 


H,-h-d 
qg= _d o 
k m (12) 
Substituting Equation 8 in Equation 12 we obtain the following expression 
m 
m, aA = (H,- h- 4d) E (13) 
The careful examination of Equations 7, 10, and 11 indicates that the 
values of A, h, and E, (at given H @’ b, ho» ds m, m,) are functions of ay 
in other words 
A= 


14 
h = fy (14) 


E= f, 


Equation 13 has only one unknown, ay the other unknowns, h, S, and 7 


have been resolved by the previous mathematical operations. Equation 13 
can now be solved by the cut and try method or by a graphical method. The 
graphical method is the much simpler and is obtained by plotting in Figure 3: 


m 
F, (a,) = m, %0 (15) 


Fy (a,) = (Hy -h- d,) E (16) 
After plotting the data of both functions between resonable limits, the 
desired value of a, will be determined as the abscissa of the intersection of 


the curves of the two functions. 
As soon as the value a, is found, the other unknowns can be determined as 


follows: The value S can be computed from Equation 5. After substituting 
a, in Equation 7, the value for A can be computed, and by using Equation 9, 
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the value of h can be found. Finally by using Equation 8, the percolation flow 
per unit length of dam can be determined, provided the coefficient of permea- 
bility is known. All above discussion applies when hy > 0. When ho = 0, the 


equations can be used in a similar manner as follows: 
Equations 2 and 12 remain unchanged, and basic Equations 3, 4, and 5 
reduce to the following: 


(3.1) 
(4.1) 


(5.1) 
Also Equation 10 can be reduced and is 


-a) a? 
(10.1) 


b 
hy = 2a, (m, +H 


where hy indicates the special case. 


Substituting in Equation 12 the value of 2 from Equation 4.1 we have 


am 
m 
1 
which can be solved by cut and try or graphical means as for the case in 
which hy, > 0. After determining the value of ay S can be found by Equation 
5.1, the value hy by Equation 10.1, and from Equation 4.1 percolation flow 
per unit length of dam can be determined from 
ak 


q (8.1) 


The following numerical example illustrates the computations where h, 


is greater than zero, showing in Table 1 an arrangement of the mathematical 
expressions which simplify the computation. A similar arrangement is 
easily set up where h, is equal to zero. 


Numerical Example 


Given: Height of the dam Ha = 100.0 feet 
Freeboard d, = 10.0 feet 
Depth of water in reservoir Hy = 90.0 feet 


Width of the crest b = 30.0 feet 
Upstream slope Ctn = m = 3.00 
Downstream slope Ctn, =m, = 2.75 


Depth of tail water = hy = 45.0 feet 


Coefficient of permeability k = 2.0 feet per year 
(from laboratory tests) 
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Required to find: 


(1) The vertical distance from the surface of tail water to the outlet 
point (M,) of the line of saturation. . 


(2) The horizontal distance from Section W-W (see Figure 1) to the 
outlet point (M,).-... 


(3) The maximum ordinate at the intersection of the line of saturation 
with Section W-W ...... 


(4) The percolation discharge through the unit length of the dam. q 


Note: The numerical values are taken for this example only and should 
not be considered as recommended values for any design. 
First, we must determine the value a,’ For this purpose, use Equation 13 


and substitute different values of ay until the left function, Fi (a,), is approxi- 


mately equal to the right function, Fy (a,), as is shown underlined in the right 


hand column of Table 1. This gives the required value for as: 
The value of a, may be found also by graphical method by plotting values 
of F, (a,) and Fy (a,) from Table 1 against the value of a, in Figure 3 and 


observing the abscissa of the intersection of the two curves. 


| 

| 
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After determining a,» we can proceed to determine the other unknowns. 


By substituting in Equation 5, we have 
S = 30.0+ 2.75 [100.0 - (7.0 + 45.0)] = 162.0 feet 


and from Equation 9 we get 


2 x 162 x 7.0 x 3.01 


+ (1.0 + 45.0)" = 72.02 feet 


h = 


by Equation 8 
g_ 
k 2.75 
q = 17.66 x 2.0 = 15.3 cubic feet per year per linear foot of 
dam or 0.31 gallon per day per linear foot of dam 


x 3.01 = 7.66 and 


From the calculated values of h and . , by using Equation 1 the trace of 


the line of saturation can be plotted for the central part of the dam between 
the points M, and M,: 


For this numerical example, the equation of the curve is 


= (72.02) - 2 (7.66) x 
= 5186.88 - 15.32 x 
The plotted data are shown in Figure 4. 


Author’s note: In the designing of an earth dam, two types of movement of 
moisture in soil can be considered. These are: (a) The movement of water 
which is dependent on pore space (in saturated material) and induced by gravi- 
ty. (b) The movement of water in unsaturated soils due to capillary forces 
or water-film movements which are dependent upon capillary action of the 
soil, and the magnitude of these forces is controlled by physical and physico- 
chemical properties of soil. 

Positive statement on the influence of the “capillary rise” on the position 
of phreatic line and upon the seepage flow through a homogenous earth dam 
may prove to be premature at this time. Selected references on these sub- 
jects are listed below. 

References: 

Physics of the Earth IX, Hydrolcgy, edited by Oscar E. Meinzer 

(a) Karl Terzaghi, Chapter IX Soil Moisture, “Soil Moisture and 
Capillary Phenomena in Soils,” page 331 

(b) Leonard D. Baver, “Retention and Movement of Soil Moisture,” 
page 364 


R. K. Schofield, “The pF of the Water Soil,” Transaction of the Third 
International Congress of Soil Science, Oxford, England, 1935, page 37 


Karl Terzaghi, “Theoretical Soil Mechanics,” May 1944, Chapter XIV, 
page 297 


Donald M. Burmister and Ek- Khoo Tan, Proceedings A.S.C.E., February 
1948, Discussion, Stability of Soil Slopes, page 221 


T. E. Buchrer and M. S. Rose, “University of Arizona College of Agricul- 
ture,” June 20, 1943, Studies in Soil Structures versus Bound Water in 
Normal and Puddled Soils 
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N. J. Denisov, “Settlement Properties of the Loessial Soils.” Government 
Publishing Office. “Soviet Science” U.S.S.R. Moscow, 1946 (Abstracted 
and translated by K. P. Karpoff and H. J. Gibbs, Bureau of Reclamation, 
Denver, Colorado, June 1951, Chapter IV, pages 7 to 9). 


STABILITY OF SLOPES 


The analysis of soil slope stability developed by N. N. Pavlovsky* is based 
on consideration of fortes exerted by the percolating water. During the per- 
colation through an earth dam, the water meets hydraulic resistance, and the 
flow will spend some part of its energy in overcoming the resistance. Sucha 
condition exists in the dam, where the line of saturation intersects the down- 
stream slope (M,); and kinetic energy, not used to overcome resistance to 


flow, will exert a hydraulic pressure on the soil particles below M,- 


Pavlovsky’s theory is that the hydraulic pressures which occur in the down- 
stream region below the line of saturation of the earth structure have a defi- 
nite effect on the stability of slopes and should be taken into consideration. 

The initial failure of the downstream slope of an earth dam or earth canal 
slope as a result of ground-water pressure is usually localized in a definite 
region, i.e., below the intersection of the line of saturation with the slope, 
and this is followed by sliding of a mass of earth material above that region. 

In this article is submitted a practical method used with good results in 
the field to compute the stability of slopes from a simplification by the author 
of N. N. Pavlovsky’s theory.° The mathematical formulas involved and two 
illustrative numerical examples are given. 


Discussion 


After determining the position of the line of saturation in an earth dam by 
the method of N. N. Pavlovsky, discussed in the first part of this article, the 
stability of the downstream slope is studied by an analysis of the forces 
which act on a particle of soil at the outlet point of the saturation line. 

Two forces act on a soil particle at the outlet point (Figure 5). 


a. The weight of the particle which is directed vertically downward 


b. The outward force of the seepage flow which has a direction tangential 
to the line of saturation at the outlet point 


4. Pavlovsky, Professor N. N., “Filtration of Water through Earth Dam,” 


Nauchno-Meliorasiony Institute, 1930, (in Russian) or thesis of K. P. 
Karpoff for B. S. in Civil Engineering, April 10, 1934 Engin. Library 
University of California. 

5. Author’s note: For simplicity of mathematical analysis a vector t wo 


(Figures 5 and 6) was plotted along the slope; whereas Pavlovsky plotted 
a vector tangent ( 2) to the phreatic line at outlet point. By this 


simplification the final equation obtained was identical to Pavlovsky’s 
equation. 


| 
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These forces can be considered as acting on a unit volume of the soil mate- 
rial. The weight of this unit volume, in submerged unit weight, will be de- 
noted by J; then, according to Terzaghi’s formula, § 


T= (s- 1) (1-n) (16a) 
where 


Ss = specific gravity of soil material and 
n = porosity of material 


The outward force of the seepage water which acts on the unit volume will 
be 


Ty J 
where J is the average hydraulic gradient for the unit volume, and yy is the 
weight of the unit volume of water. 
If the resultant force R (Figure 6) of the two above-mentioned forces is 


projected on the plane of the slope and the plane normal to that slope respec- 
tively, the two components R, and R, will be obtained. Then the equation of 


equilibrium, at which the particle of soil will overcome the force of friction 
and will be ready to slip down the slope, may be written: 


R,=R, tan (17) 


where ¢ is the ai of internal friction of the soil material. 
By projecting the force T'on the plane normal to the slope (Figure 6), we 
obtain the vector Ry and 


T'Cos ( 18) 


Force ri and the component of ¥ ‘(Figure 6) on the plane of the slope, ina 


state of equilibrium, must be equal to the force of friction, Equation 17, and 
may be expressed: 


R,tan@ = J+ Sina (19) 


After substituting in Equation 19 the value R, from Equation 18 and solving 
for the limit of J, we have 


(tana - tan@) Cos a ( 20) 


y" 


w 


Jim = T"(tan? - tana ) Cos ( 21) 


6. Theory of Terzaghi--“Erdbaumechanik auf Bodenphysikalisher Grundlage,” 
1925, pp. 128-131. 
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According to Pavlovsky’s theory, the value of the hydraulic gradient (J) at 
the outlet point M, (Figure 7) of the line of saturation is equal to the slope of 


the dam; that is, it coincides with the slope of a structure, or in this case 


y 
must represent a true gradient (J — 
0 


terial which lies on the slope to be stable, J 
or 


tana = ). In order for a soi) ma- 


must be greater than J 


lim true’ 


T(tan ¢ - tana ) Cosa> tana (22) 
Now introducing the safety factor “c” to Equation 22 we obtain 
Y"(tang - tana ) Cosa@ =c tana (23) 


This expression represents the largest possible value of hydraulic gradient 
for stability at the point under consideration. 

The stability of an earth slope can be improved by placing a blanket of 
selected granular material on top of the slope, and in this case, Equation 23 
is expressed as follows: 


(Y¥"+as') (tan@ - tana) Cosa =c tane (24) 


that is, the stability of slope is increased by Y"+as'or1+ as' times where 
a = the thickness of the layer of granular material 
(normal to ihe slope) 


s'= bulk specific gravity of the granular material 


Instead of computing tana@ by using Equation 24, it is possible to use a graph 
(Figure 8) which is plotted on the basis of Equation 24 where c and ¥ “are 
considered equal to unity and the thickness of protective blanket from granu- 
lar material a = 0. Obtaining the value p§ from the graph (Figure 8), we can 
determine the true value for the downstream slope of the structure by using 
Equation 25 


D c 
Slope = (25) 


For the practical application of the formula which is developed above, it is 
only necessary to have the following information about the soil: 


Angle of internal friction (¢) of soil in a saturated condition 
Specific gravity of soil 

Average dry density 

Bulk specific gravity of the protective granular material 


Numerical Example 


Given: (1) Material, clay type 
(2) Specific gravity - 2.71 
(3) Angle of internal friction in a saturated 
condition tan ¢ = 0.424, ¢ = 23° 
(4) Average dry density - 106.0 pounds per 
cubic foot 
(5) Solid density = 62.4 x 2.71 = 169.1 pounds 
per cubic foot 
7. 1 =(tan@ - tana) Cosa 
D 


| 
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Then porosity, n = 190.1 we om . 37.3 percent, and weight per unit 


volume of soil at the hydraulic pressure 


¥'=(s- 1) (1-n) = (2.71 - 1) (1 - 0.377) = 1.07 grams 
per cubic centimeter and 


T"= 1.07 


From the graph (Figure 8), for @ = 23°, the D = 4.75 
Using 1.1 for the safety factor and introducing it in Equation 25, we obtain 


Slope = | :1 =4.88:1 


Hence, the slope of the downstream face of the structure for the sake of sta- 
bility should be about 5:1. 

Let us see what will happen if a layer of selected granular material 0.458 
meter is placed on this slope. Bulk specific gravity of the granular material 
(s') = 2.60, then substituting in Equation 25 


1.1 
1.07 + 0.458 x 2.60 


and the allowable slope of the structure with a protective blanket is about 
2-1/4:1. This example shows clearly the utility of loading the lower portion 
of the downstream part of the earth dam with a pervious material. 


Note: The safety factor c can vary from 1.1 and 1.3. The selection of a 
factor during the computation of slope depends upon the physical properties 
of soil material, field conditions, and the judgment of the designer. 

It is believed that the above method, due to its simplicity and ease in 
mathematical approach, is of value for designers of earth hydraulic struc- 
tures. It is particularly valuable to the field engineer where a quick decision 
must be made for determining the stability of a structure. 


Slope = [4.75 
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Line of saturation 


FIGURE 5 
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c tana 
ton a) Cos a 


=(tan §- tona)Cos ¢ 
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(tan 
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function of % ond this relationship is 


Above expression shows that Cot ¢ or Dis 
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